The final stage of gravitationally collapsed thick matter layers 
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In the presence of a minimal length physical objects cannot collapse to an infinite 
density, singular, matter point. In this note we consider the possible final stage of 
the gravitational collapse of "thick" matter layers. The energy momentum tensor 
we choose to model these shell-like objects is a proper modification of the source 
for "non-commutative geometry inspired", regular black holes. By using higher mo- 
menta of Gaussian distribution to localize matter at finite distance from the origin, 
we obtain new solutions of the Einstein's equation which smoothly interpolates be- 
tween Minkowski geometry near the center of the shell and Schwarzschild spacetime 
far away from the matter layer. The metric is curvature singularity free. Black hole 
type solutions exist only for "heavy" shells, i.e. M > M e , where M e is the mass of 
the extremal configuration. We determine the Hawking temperature and a modified 
Area Law taking into account the extended nature of the source. 

PACS numbers: 04.40.-b, 04.70.Dy, 04.50.Kd 



o 



5-1 . 

bo 



> 
(N 
m 



Keywords: 



INTRODUCTION 



Relativistic, self-gravitating matter shells have been throughly investigated in different 



CO ■ sectors of theoretical physics: "... from cosmic inflation to hadronic bags" Remark- 

able applications of gravitational shell models can be found in the framework of inflationary 
cosmology, where, both the "birth" and the evolution of vacuum bubbles can be effectively 
described in terms of the dynamics of the boundary surface engulfing a false vacuum do- 
main Matter shells in general relativity are modeled as "zero-thickness" membranes 
endowed with some characteristic tension determined by the underlying classical, or quan- 
tum, physics. Neglecting the real width of the mass-energy distribution affects Einstein 
equations by introducing a surface of discontinuity in the background spacetime. This ap- 
proximation allows to encode all the dynamics in the matching condition between the inner 
and outer geometries Furthermore, contracting matter shells provide useful analytic 
toy-models of collapsing massive bodies leading to black hole formation. In the spherically 
symmetric case the only dynamical degree of freedom is given by the shell radius and the 
system can be quantized according with the standard principle of quantum mechanics. In 
this framework, self-gravitating quantum shells open a window over the still "murky" quan- 
tum features of evaporating mini black holes 



In this paper we are going to investigate the static, final stage, of collapsed matter shell in 
the presence of a fundamental minimal length forbidding the shell to contract into a singular 
matter point. The emergence of a minimal length, as a new fundamental constant of Nature 
on the same ground as c and h, is a general feature of different appraches to quantum gravity 
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nHi3 . 

In a series of papers we showed as the very concept of "point-particle" is meaningless if 



there exists a lower bound to physically measurable lengths [14H18I]. This basic notion can 



be encoded into the Einstein equations through a proper choice of the energy momentum 



tensor |19l-l29l|. The most remarkable outcomes of this procedure are the disappearance of 
curvature singularities in the solutions of the Einstein equations j30|, a regular behavior of 
the Hawking tem per ature which allows to determine the physical character of the evapora- 
tion remnant [31-37j], and a different form of the relation between entropy and area of the 
event horizon reproducing the celebrated Area Law for large, semi-classical, black holes. 
The paper is organized as follows. In section [Til we will derive new solutions of Einstein 
equations describing, spherically symmetric, static, self-gravitating, thick matter layers, in 
the presence of a fundamental minimal length. These type of objects are modeled by means 
of a proper modification of the source for "non-commutative geometry inspired", regular 
black holes. We replace the Gaussian profile of mass-energy, peaked around the origin, with 
higher moments of Gaussian distribution with maxima shifted at finite distance from the 
origin. The finite width of the distribution is determined by the minimal length. In the case 
of lump-type objects, the minimal length measures the spread of mass-energy around the 
origin and removes the curvature singularity. For finite width matter layers, with energy 
density vanishing at short distance, not only the central curvature singularity is removed, 
but the extrinsic curvature discontinuity between "inner" and "outer" geometry is cured, as 
well. Spacetime geometry is continuous and differentiable everywhere and smoothly inter- 
polates between Minkowski metric near the center, and Schwarzschild spacetime far away 
from the matter layer. 

As in the cases previously discussed, black hole type solutions exist only for "heavy" shells, 
i.e. M > M e , where M e is the mass of the extremal configuration. "Light" matter layers 
with M < M e will settle down in a smooth solitonic type configuration with no horizons or 
curvature singularity. 

In section [TTT1 we study the thermodynamic properties of black hole solutions and determine 
both the Hawking Temperature, T#, and the relation between entropy and area of the hori- 
zon. We recover the celebrated Area Law in the limit of large, semi-classical, black holes. 
On a general ground, we find the leading term is one fourth of the area but in units of an 
"effective" gravitational coupling constant, Gat (?"+), depending from the radius of black 
hole. For "large" r+, the Newton Constant is recovered. Finally, in section [IV] we will draw 
the conclusions. 



II. THICK SHELLS 

In a previous series of papers we solved the Einstein equations including the effects of a 
minimal length in a proper energy-momentum tensor 

R^u — -g^u R = 8tt G N T^ U , (l) 

where = diag(p ,p r ,p± ,p±). The energy density was chosen to be a minimal width 
Gaussian distribution 



3 



representing a "blob-like" object, centered around the origin, with a characteristic extension 
given by I oc \J~Q. The radial pressure p r is fixed by the equation of state p r = —p reproducing 
the deSitter "vacuum" equation of state at short distance. This is a key feature to build up 
a regular, stable configuration, where the negative pressure balances the gravitational pull. 
In other words, the singularity theorem is evaded by a violation of null energy condition 
triggered by the short distance vacuum fluctuations. 

Finally, the tangential pressure p± is obtained in terms of p by the divergence free condition 
T^ u = 0. 

As the source is static and spherically symmetric, the line element can be cast in the form 

ds z = -f( r ) dt * + ^ + r 2 n 2 , (3) 



with 



f(r) = 1 - 2GjV r m(r) - (4) 



The cumulative mass distribution m(r) is given by 

m {r) = 4tt / dr\r'f p (r'). (5) 
Jo 

We notice that the parameter M corresponds to the total mass energy of the system, namely 

M = lim m(r). (6) 

r— >oo 

The above profile cures the usual Dirac delta (singular) distribution associated to a point- 
particle, and leads to a family of regular black hole solutions [l9l-28 



We are going now to adapt the above approach to an object which is shell-like, rather than 
point-like. For this reason, we recall that the Gaussian distribution admits moments of all 
orders. The regular black holes we found from Gaussian distribution are just gravitational 
objects corresponding to the zero-th moment of the distribution. Here we want to explore the 
nature of the solutions emerging from higher moments of the Gaussian, namely normalized 
matter distributions of the kind 



r 2fc e -r 2 /46> 



p k (r)-M 4fc+2?r fc+ 3 /2 r ( ^ + 3/2 ) (7) 

where k = , 1 , 2 . . . is a natural number. We are considering the case of even moments 
only. This choice is motivated by the request of having a rapidly vanishing energy profile at 
the origin, corresponding to the case of "empty" shells. On the contrary odd moments may 
admit a linearly vanishing energy profile, which gives rise to a gravitational object which 
does not properly describe a matter shell. We stress that the choice of even moments does 
not modify any of the following conclusions for both the geometric and the thermodynamic 
properties of the solutions. 

For k = the function p k turns into the Gaussian distribution, centered around the 
origin, while for k > 1 the matter distribution is more and more diluted near the origin, 
being peaked at tm = 2yk9. As a result the density function (jTj) describes a whole family 
of "mass-degenerate" shells, with the same M, but concentrated at a distance given by r« 
(see Fig. [Q. Eq. (j7|) discloses further properties of the minimal length, which assumes a 




FIG. 1: Plot of the radial density profile p(r) as a function of r in VO-units for k = 1 (cyan), 
k = 2 (blue), k = 7 (red), k = 8 (pink), k = 20 (green), k = 21 (yellow). The function p(r) has 
been normalized so that its integration over a spherical volume gives 1. 



new intriguing meaning. In the case of point-particle, y9 represents the spread of the object 
around the origin. For matter shells, we see \/6 relates to the shell thickness or, in other 
words, a measure of the intrinsic fuzziness of the layer. Just as a particle cannot be exactly 
localized at a single point, we cannot have zero-width layers, as well. Thus, as the matter 
distribution is smooth everywhere, there is discontinuity in the extrinsic curvature between 
the "inner" and outer "geometry". No matching condition is required and we can look for 
a single, smooth, metric inside, across and outside the matter layer. 

We can define the distance between two shells corresponding to different moments as the 
distance between the peaks: 



2V0 



Vk + l-Vk . 



(8) 



We see that for higher moments Ar^/ vanishes as we are at length scale much larger than 
\/9 and the relative distance cannot be resolved anymore. 

A stable solution of the Einstein equations can be obtained by sourcing the gravitational field 
by the energy momentum tensor of an anisotropic fluid: the choice p r = —pt for the matter 
equation of state allows to have goo = ~9rr ■ Furthermore, the hydrodynamic equilibrium 
equation will give pt in terms of pk, while the metric itself results to be independent from 
p t . The solution of the Einstein equations reads in geometric units, c = 1, — 1: 



ds 7 



2m ( r 



r 



dt 2 + 1 - 



2m ( r 



r 



dr A + r W 



7 (A; + 3/2 ;r 2 /46 



m(r)=An / dr' r' 2 p (r' ) = M ■ 

K J Jo T(k + 3/2 



(9) 
(10) 



where 7 ( & + 3/2 ; r 2 /49 ) is the lower, incomplete, Euler Gamma Function. By taking into 
account the asymptotic form of 7 ( k + 3/2 ;r 2 /A6 ) for small argument, one finds that the 
metric is essentially flat near the origin, the contrary to what happens for the case k = 
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FIG. 2: Plot o/goo( r ) as a function of r in \fd- units for fixed value of k = 1 and different values 
of the mass M , i.e. M = 1 (blue), M = M e (magenta) and M = 5 (brown). Here M e ~ 2.36976 
is the extremal mass for the case k = 1 . 



in which a regular deSitter core forms. Indeed one finds 



f(r) 



1 



2k + 3 



4M 



_rT(A; + 3/2 



/ r 2 \ k+3/2 



<is 2 = -(l-0(r 2M )) dt 2 +(l-0(r 2M )y 1 <ir 2 + rW . (12) 

The asymptotic behavior f[T^]) can be seen as a generalization of the "Gauss Theorem" : 
inside an empty, classical, thin shell of matter the Newtonian gravitational field is zero, i.e. 
spacetime is flat. Introducing a finite width density, smoothly decreasing both toward the 
origin and space-like infinity, causes small deviations from perfect flatness. 
Going back to (J9j), the presence of event horizons is read from the zeroes of g~} (see Fig. [2]). 
Due to the non-trivial structure of the metric, it is more convenient to look at the "horizon 
equation" in the form 

UsU{TB) , u{rB)s ^ -m^— . (13) 

The problem of finding the horizons in the metric ([9]) is mapped in the equivalent problem 
of determining the turning-points for the motion of a "test-particle" of energy M subject 
to the "potential" U (r H ). The intersection (s) between a line M = const, and the plot of 
the function U(th), in the plane M-r^, represent the allowed radii of inner/outer event 
horizons (see Fig. [3]). The asymptotic behavior of U ( r# ) is obtained from the corresponding 
approximate forms of 7: 

U{r H )^Y' r H»^ (14) 

/ 2 \ 

u(r H )^Ver(k + 5/2) r^j , r H «V6. (15) 

Thus, U ( th ) is a convex function with a single minimum determined by the condition 



dU „ r 



2fc+3 



-> o 2fc+ L fc+ 3/ 2 =7(^ + 3/2; rl/AO ) . (16) 



dr H 2 2k + 2 6 k + 3 / 2 
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FIG. 3: Plot of U(th) cls a function of ru in VO-units for k = (blue), k = 1 (magenta) ant 
k = 4 (brown). 



The radius r e represents the size of an extremal configuration with a couple of degener- 
ate horizons: r_ = r + = r e . Once r e is numerically determined from ffTB]) . one gets the 
corresponding mass M e from the potential: 



M P = U( r P 



(17) 



An order of magnitude estimate for r e and M e can be obtained by keeping only the leading 
theta-term: 

r P ocVd — > M e oc - — M PL (18) 



L 



pi. 



where Lpi is the Planck length, i.e., Lpi = Mp X = \JGn- This estimate suggests that 
(near/) extremal configurations are close to a full quantum gravity regime. As such, they are 
appropriate candidates to describe the end-point of the Hawking evaporation process where 
the semi-classical description breaks down. 
In summary: 

• for M > M e we find a geometry with non-coincident inner and outer horizons r_ < r + . 
It is worth to remark that there is no curvature singularity in r = 0. Spacetime is flat 
near the origin. 

• For M = M e we have an extremal configuration with a pair of degenerate horizons, 
r_ = r + = r e . 

• For M < M e there are neither horizons nor curvature singularities. The metric is 
smooth and regular everywhere. The shell is too light and diluted to produce any 
relevant alteration of the spacetime fabric. After collapse it will settle into a sort of 
solitonic object with no horizon of curvature singularity. 

Massive layers will produce horizons, but no curvature singularities. This is a crucial point 
which marks a departure with respect to all the existing literature. Our approach must 
not be confused with previous contributions in which generalized matter shells containing 
polytropic and Chaplygin gas cannot ultimately resolve the emergence of singularities [38] . 
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We can estimate the size of these objects by solving iteratively the equation ( ITS]) and 
truncating the procedure at the first order in the expansion parameter exp ( —M 2 /6): 



■>'+ 



2M 



1 _ (M 2 /^+V 2 e _ MVe 



r(fc + 3/2) 

This quantity can be compared with the shell mean radius 

T(k + 2) 



47T 

M 



/■oo 

/ drr 2 r p(r) = 2\f§ 
Jo 



T(k + 3/2) ' 

We notice that for k » 1 the leading contribution to the ratio r+/(r ) is given by 

r, M 1 



(19) 



(20) 



(21) 



Thus, for an assigned M the ratio decreases with k. The horizon radius is determined by 
the total mass energy M and is weakly fc-dependent. On the contrary, (r) oc y/k + 1 a/# 
and grows with fe, which means that for higher moments the horizon is surrounded by a 
cloud of matter. 



III. THERMODYNAMICS 

Massive shells will collapse into black holes described by the line element (jSD,®- This is 
not the end of story as these objects are semi-clasically unstable under Hawking emission. 
We are now ready to study the thermodynamic properties of these solutions starting, from 
the Hawking temperature: 



H 47rr_i 



1 



r 2fc+3 e -r 2 + /Ae 
2 2fc+2#fc+3/2 ~ ( I + 3/2 j r 2 / A Q ) 



(22) 



It can be easily verified that T# is vanishing for the extremal configuration Tjj ( r + = r e ) = 0. 
For k = we obtain the temperature of the black hole as in [2lJ. The behavior of the 
temperature can be found in Fig. HI We notice that the regularity of the manifold for any k 
leads to a cooling down of the horizon in the terminal phase of the Hawking process. At r + = 
r max the presence of a maximum temperature corresponds to an infinite discontinuity the 
heat capacity which is usually interpreted as the signal of a "change of state" for the system. 
For r > r max the black hole is thermodynamically unstable and increases its temperature by 
radiating away its own mass. After crossing r max , i.e., for r + < r max the black hole enters 
a stability phase asymptotically ending into a degenerate ( zero-temperature ) extremal 
configuration. By increasing k we just lower down the maximum temperature. As a result 
the solution is unaffected by any relevant quantum back reaction as already proved for the 
case k = in (2ll |. 

A further important consequence of Th < T max is that ( in extra-dimensional models ) the 
black hole mainly radiates on the brane "j43[ , it never becomes hot enough to warm-up the 
bulk in a significant way. 

The area-entropy law can be recovered from the relation 

dM = T H dS (23) 
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FIG. 4: Plot o/T(r+) as a function of r + in y/9 -units for k = (blue), k = 1 (magenta) and k = 4 
(brown), k = 15 (green). The dashed line represent the classical case, i.e. 9 = 0. 



TABLE I: Some values of parameters of the matter shells 







M e 


T 

- 1 max 


''max 


k = 


3.0224 


1.9041 


0.014937 


4.76421 


fc = 1 


4.0431 


2.3698 


0.012783 


5.72632 


k = 4 


5.9269 


3.2647 


0.009960 


7.56069 


fc= 15 


9.7347 


5.1245 


0.006799 


11.3419 



where dS is the horizon entropy variation triggered by a variation dM in the total mass 
energy M. In order to translate the first law of black hole thermodynamics (1231) into a 
relation involving the area of the event horizon, we need to write mass energy variation as 

dU 

dM=— dr + (24) 
or + 

and to take into account that the minimum of U ( r+ ) is the mass of the extremal black 
hole. Thus, when integrating ( 124"|) the lower integration limit is the radius of the extremal 
configuration 

pre I QU pre 

S= dx—— = 2nT(k + 3/2) dr - f - ■ — -— — , (25) 



T H dx v Jr + l(k + 3/2; r 2 /46 , 

where we have inserted (j22p and (j24p into (1251) . By performing the integral one gets 



S= nT(k + 3/2 



1 9 

rl rz 



7 ( k + 3/2 ; r 2 + /40 ) 7 ( A; + 3/2 ; r e 2 /40 ; 



7rr(A; + 3/2) / drr^. (26) 
The first term can be written in terms of the area of the event horizon as 

S = 4G^ T{k + 3/2) (7(^ + 3/2%^) "7(^ + 3/2;r e 74£)) + --- (2?) 
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and represents the "Area Law" in our case. We have re-inserted the Newton constant into 
( 127)) for reasons to become clear in a while. The standard form, which is one fourth of the 
area, is recovered in the large black hole limit, r# >> y/9. 

Once ()27j) is written in natural units, we can define an effective Newton constant as 

7 (A; + 3/2 ;rl/46) 
G N ^G N (r + ) = G N T(k + 3/2) (28) 

and introduce a "modified" Area Law as 



S(r H ) = -^ (29) 
AG N ( r H ) 

The interesting feature of Gn(th) is to be "asymptotically free" in the sense that it is 
smaller than Gjy, which represents the asymptotic value of the gravitational coupling for 
large black holes only j48[. 



Finally, the second term in ( )26|) . gives exponentially small corrections to the leading area 
term. It can be analytically computed for k » 1, but it does not introduce any relevant 
new effect. 



IV. CONCLUSIONS 

In this paper we derived new static, spherically symmetric regular solutions of Einstein 
equations, describing the final state of collapsing matter shells in the presence of an effective 
minimal length. This derivation is the result of a long path started with the quest of quantum 
gravity regularized black hole solutions (30[. We showed that our previous discoveries of 
noncommutative geometry inspired solutions correspond to the simplest case within a larger 
class of regular gravitational objects. This class, the family of gravitational shells here 
presented, has peculiar properties which descend from the common key feature of the absence 
of curvature singularities. Such properties are: the possibility of horizon extremisation even 
for the neutral, non-rotating case, the existence of a minimum amount of mass energy below 
which horizons do not form, the occurrence of a phase transition from a thermodynamically 
unstable classical phase to a thermodynamically stable cooling down in the final stage of the 
horizon evaporation. As for the case of regular black holes, these matter shells correspond to 
final configurations of dynamical processes of gravitational collapse in which matter cannot 
be compressed below a fundamental length scale. The above properties are also common to 
other models of quantum geometry, that have been derived by means of different formalisms. 
Specifically the shell profile of the energy density resembles what one has in the case of loop 



quantum black holes, with consequent thermodynamic similarities [44J-I471 . 

The study of these objects is far from being complete. We believe that these new solutions 
can have repercussions in a variety of fields. Here we just mention that these re gula r shells 
could affect the stability of the deSitter space, at least during inflationary epochs [27| as well 
as they could lead to new insights in the physics of nuclear matter via the gauge/gravity 
duality paradigm [33]. 
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